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The reahzation of the spin-Hall effect in quantum wells has led to a plethora of studies regarding 
the properties of the edge states of a 2D topological insulator. These edge states constitute a class of 
one-dimensional liquids, called the helical liquid, where an electron's spin quantization axis is tied to 
its momentum. In contrast to one dimensional conductors, magnetic impurities — below the Kondo 
temperature — cannot block transport and one expects the current to circumvent the impurity. To 
study this phenomenon, we consider the single impurity Anderson model embedded into an edge 
of a Kane-Mele ribbon with up to 512 x 80 sites and use the numerically exact continuous time 
QMC method to study the Kondo effect. We present results on the temperature dependence of the 
spectral properties of the impurity and the bulk system that show the behaviour of the system in 
the various regimes of the Anderson model. A view complementary to the single particle spectral 
functions can be obtained using the spatial behaviour of the spin-spin correlation functions. Here 
we show the characteristic, algebraic decay in the edge channel near the impurity. 

PACS numbers: 71.10.Fd, 05.10.Ln, 05.70.Ln 



I. INTRODUCTION 

Numerical studies of variants of the Kane-Mele modeP 
have recently been pursued with increasing intereslP^^ 
since it can be used as a theoretical framework to study 
correlation effects in quantum spin-Hall insulators, or 
two-dimensional (2D) topological insulators (TIs)^. A 
characteristic feature is the formation of metallic edge 
states at the boundary of the system, which are robust to 
external perturbations provided that time reversal sym- 
metry is not broken.l^ They form a helical liquid, such 
that the electrons' spin is tied to its direction of motiorP. 

A particularly interesting perturbation of the helical 
edge state is the introduction of magnetic impurities in- 
teracting with the edge — a problem usually modeled by 
an S' = i local spin that is coupled to the helical liq- 
uid. Due to the one dimensional(lD) nature of the edge 
this problem has been studied extensively with bosoniza- 
tion techniques with variable Luttinger liquid parame- 
ter accounting for correlation effects in the helical edge 
state^"^^''. The Kondo effect in three dimensional TIs has 
been studied in Refs.'^^Hlll. Returning to the 2D case and 
in the weak coupling regime with respect to electronic 
correlations on the edge, the formation of the Kondo sin- 
glet will effectively remove sites, thereby redefining the 
topology of the slab and the flow of the edge state. The 
aim of this paper is to study the temperature dependence 
of this effect by computing, among other quantities, the 
site dependent density of states. To do so, we will set 
out to model the magnetic impurity using the single im- 
purity Anderson model^^ which accounts for a single lo- 
calized energy level that hybridizes with the states of the 
TI and has an on-site Coulomb repulsion, while sites in 
the bulk are assumed to be noninteracting. For non- 
vanishing Hubbard-interaction this model enables us to 
trace the progression from the high-temperature regime 
over the development of the local moment towards the 



formation of the Kondo singlet. We achieve this with 
the numerically exact interaction expansion continuous 
time Quantum Monte Carlo (CT-INT) algorithm intro- 
duced by Rubtsov et al.'iSttSSI which is particularly suit- 
able for the study of impurity problems, since nonin- 
teracting bath sites (the TI) do not count towards the 
computational complexity of the algorithm and can be 
integrated out. To access the single particle properties of 
the bath we calculate the self-energy on the impurity and 
then calculate the bath Green's functions using Dyson's 
equation. This allows us to exhibit the deflection of the 
edge state at the impurity by looking at the spectral sig- 
natures arising in the bulk spectral functions due to the 
emerging Kondo effect on the impurity. Towards the end 
we progress from the one-particle spectral functions to a 
two-particle quantity, the spatially resolved, equal-time, 
spin spin correlation function. Using this quantity, we 
can inquire the spatial extent of the Kondo screening 
cloud in the bulk and along the edge. 



II. THE MODEL 

The Kane-Mele model was first proposed as a candi- 
date for a possible quantum spin-Hall effect in graphene 
and its Hamiltonian is given by 



Hkm — Ht + Hx , 

with 

i7t = -i V at (b^ + b^^ - 
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FIG. 1. (Color online) The ribbon is periodic along the r - 
direction and at site r = 0, n — the orange box denotes the 
impurity orbital which couples with a matrix element V. Sub- 
lattice A is denoted by open circles and sublattice B by filled 
circles. Oi and 02 denote the unit vectors of the honeycomb 
lattice. 



The effect of electron electron interactions in graphene 
was studied in Ref.l^i^tSll, Although it turned out that in 
graphene the spin-orbit coupHng is too small to observe 
the QSH state, the model can still be used as an effective 
Hamiltonian for this topological state of matter, aj^ and 
denote fermionic operators acting on the respective 
sublattice and, to account for edge states, we consider 
this model on a slab geometry. The boundary conditions 
are periodic in r-direction - the corresponding number of 
sites is denoted by ~ and open in n-direction with a 
length of Ny sites. A is the strength of the spin-orbit in- 
teraction and the hopping t is set to t — 1 for everything 
that follows. This model is related to the spinless Hal- 
dane model that shows a quantum Hall effect but breaks 
time reversal invariance (TRlj23]. The Kane-Mele model 
can be understood as two copies of the Haldane model 
while preserving time-reversal symmetry and exhibiting 
a quantum spin-Hall effect. Into this bath system we em- 
bed an impurity at an edge. The impurity's Hamiltonian 
i?inip is given by 

^^Imp = Ho + Hu (2) 

with 

Hu^U{n^^-\){4-\). 

Here, ed denotes the energy of the dot, U the Hubbard in- 
teraction and V the hybridization between the first bath 
site and the impurity. denotes fermionic operators 
acting on the impurity. We have chosen a symmetric 
representation of the Hubbard interaction that sets the 
chemical potential to zero for the half-filled case. We 
note that the impurity Hamiltonian obeys time reversal 
symmetry together with the bath. 
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FIG. 2. Spectrum of the Kane-Mele model. Here we have 
used — 512, A^^ = 80 and A = 0.1, which constitute our 
"canonical values" in the following. Visible are the different 
bands due to the "orbitals" in n-direction, as well as the fa- 
mous edge states crossing at the fermi energy. 

III. SUMMARY OF BATH PROPERTIES 

As already mentioned, the bath model Hkm exhibits 
the so-called edge states which, as the name implies, are 
localized at the edges. Since we attach the impurity to 
a site belonging to an edge we revisit some properties of 
the bath. We refer to anything outside of the impurity as 
bath and everything in the bath that is not dominated by 
the edge state as bulk. The edge states correspond to the 
states in the energy spectrum of |Fig. 2| that cross at the 
Fermi energy, e{kx) = 0, and enable gapless electronic 
excitations at the edge. These edge states constitute a 
helical liquid where the spin of an electron is coupled to 
the direction of propagation, hence an interaction hipping 
the spin reverses its momentum. Since we will argue quite 
a bit with the help of the spectral functions we point out 
the general structure of A„(r, here. It is 

Anir,Lo) ^ A"Jlj) + Bn{r,uj,V) + Cnir,u,^iu)) (3) 

with the impurity independent background A'^(uj), a 
term i3„(r, w,]/) that depends on the hybridization 
V between lattice and impurity and the contribution 
C„(r, w, S(a;)) due to the self-energy T,{uj) of the impu- 
rity. In |Fig. 3| we show a site-resolved view onto the 
spectral functions (w) of the bath. The bulge that is 
visible in the outermost (n = 0) spectral functions is the 
edge state, which has its spectral weight centered around 
a; = 0. Since we consider the system without an impu- 
rity wc have translation invariance along the r - direction. 
For comparison we show in |Fig. 4| a cut along r = of the 
same spectral functions. Further into the bulk the gap of 
the insulator appears. Also we see the odd-even pattern 
close to the edge. The n = 1 spectral function shows a 
gap, whereas the n = 2 function shows some remains of 
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the exponentially decaying edge state. 




FIG. 3. (Color online) This plot shows a part of the real- 
space lattice at A = 0.1 where to every lattice point we have 
attached the spectral function A'^{uj). For a given value of r, 
the black arrows denote the path that is taken through the 
lattice if n is increased. Since the impurity is not yet added we 
have translation invariance along the r direction. The bulge 
of the edge state in the n = spectral functions is clearly 
visible. 
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FIG. 4. (Color online) Spectral functions of Fig. 3 for com- 
parison in a more traditional 2D plot. Both plots are from a 
system with A^^ = 512 and Ny = 80 sites at A = 0.1. The 
fine wiggles in the spectral functions, as e.g. in n = or 
n — 1, are artifacts of the finite system size. The edge state 
at n = is clearly visible which almost immediately decays 
farther in the bulk. Having used a finite rj = 4Aa; ~ 0.03 we 
have some broadening which prevents the gap of the spectral 
function from vanishing in the bulk. Alj is the resolution of 
the energy u. 



IV. AN UNCORRELATED IMPURITY 

Adding to the bath system the uncorrelatcd impurity 
given by Hq at site r = and n — the spectral prop- 
erties of the system change around the impurity since 
now the B - term in |Eq. (3) contributes to the spectral 
fimctions due to the hybridization V. From Fig. 5a and 
|Fig. 5b| it is visible that right at the impurity the spectral 
weight of the edge state is drastically reduced. In con- 
trast to the changes due to correlations in Fig. 12a we see 
that locally this potential poses quite a strong perturba- 
tion to A{uj). Since time reversal symmetry is present, 
single particle backward scattering is prohibited as this 
would amount to flipping the orientation of the spin. As 
a consequence, the edge state has to circumvent the po- 
tential impurity by deflecting into the bulk. Thus the 
missing spectral weight reappears at sites further into 
the bulk. Fig. 5c and [Fig. 5d| show the change in the 
spectral function, Bn{r,uj,V), due to the hybridization. 
The edge state is deformed around the impurity since it 
acts as a pure potential scatterer for the edge state. This 
is consistent with the spectral function of the impurity 
which is just a lorentzian around a; = similar to the 
"correlated" spectral function at /3 0.1 in |Fig. 7a| Al- 
though the edge state is protected by symmetry against 
potential scattering, the effect of the impurity is that it 
acts as a trap for electrons from the bulk which can then 
in turn interact with the electrons of the edge stat^^^l. An 
interpretation of the resulting new path of the edge chan- 
nel is that the site to which the impurity is connected is 
effectively removed from the lattice. In that sense the 
deformation of the edge state can be understood as a 
rerouting along the changed edge of the system. 



V. A CORRELATED IMPURITY 

Now we add the Hubbard interaction Hjj to the impu- 
rity which leads us to consider the single impurity Ander- 
son model (SIAM) with a bath given by the Kane-Mele 
model, 

H ~ i?Imp + Hkm- (4) 
Integrating out the bath electrons we obtain the action 

S = -Y,jdT j dT'dl{T)GaAr-r')dAT') 







(5) 



+ U I dr{n''^{T)~\){n'l{T)-\) 



which is perfectly suitable for an implementation of a 
numerically exact CTQMC method that expands in the 
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FIG. 5. (Color online) (a) shows the complete spectral functions now with a non-interacting impurity located at r = and 
n = Q. The original edge state is deformed into the bulk around the impurity. The missing weight in the n = spectral 
functions shows up in the n = 1,2 spectral functions, (b) is a cut of y4„(r, tj) along r = that shows the rearrangement of 
spectral weight from the edge into the bulk, (c) shows Bn{r, uj,V), the effect on the spectral functions attributable to the 
hybridization V of the impurity with the edge state. We see large negative contributions to the edge state at the impurity and 
positive contributions farther into the bulk. Bnir = Q,UJ, V) is shown in (d). Obvious is the strong reduction at the site below 
the impurity (r = 0, n = 0) which is shifted to the sites around the impurity. 



interaction strength U . The computational effort is re- 
duced if one uses the property that our model is time 
reversal invariant which leads to a spin-diagonal Gd.d ^ 
detailed in |Sec. B| for a general time reversal symmet- 
ric Hamiltonian. To at best study the temperature de- 
pendence of the above impurity problem, we will com- 



pute several quantities 
in 



Fig. 



The double occupancy, (n^njj) 



6a will allow us to track the formation of the lo- 
cal moment and its screening. The same information is 
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FIG. 6. (Color online) Overview of, (a), the double occupancy 
and in (b) , Xzz , as a function of temperature T — ^. The plots 
of the double occupancy highlight the different regimes. All 
functions start out at the uncorrelated value {nfn_i) — 0.25 
and then fall towards some dip at intermediate temperatures. 
This dip roughly coincides with the local moment regime. 
After that, the double occupancy increases slightly up to a 
satiiration which is a signature of the Kondo regime. 






FIG. 7. (Color online) The spectral functions of the dot for (a) 
U — 2 and (b) (7 = 4. At /3 = 0.1 we see the high-temperature 
regime with weight centered around zero. Lowering the tem- 
perature we cross over into the local moment regime with 
two clearly seperated Hubbard bands at oj « Lowering 
the temperature further we see the emergence of the Kondo 
resonance at = 0. 



essentially contained in the local spin susceptibility 



Xzz = / dT{S^{T)S^) 



measured on the impurity and plotted in |Fig. 6bl Due to 
TRI it is sufficient to consider only Xzz since the other 
components are degenerate. 

To study the deflection of the edge state as well as the 
formation of the Kondo resonance, we also compute the 
site dependent single particle density spectral function 



;Im(G„(r, w + i?7)) 



(6) 



for a small t]. To obtain the spectral functions of the 
bath lattice we analytically continue the impurity spec- 
tral function to real frequencies and use Dyson's equa- 
tion to access the lattice spectral functions as detailed 
in appendix [A] Additiona l ly we g et access to the term 
C„(r, w, E(aj)) of Eq. (3) Fig. 8 shows the self-energy 



FIG. 8. (Color online) Real and imaginary part of the self- 
energy E(tj) in different temperature regimes for U — 2. 
Im(E(cj)) shows a simple lorentzian shape for /3 = 2 in the 
local moment regime. Crossing over to the Kondo regime we 
see the development of a two peak structure. Since 'S{z) is a 
holomorphic function its real and imaginary part are linked 
via the Kramers-Kronig relations. 



different temperature regimes. Since ^{z) is a holomor- 
phic function its real and imaginary part are linked via 
the Kramers-Kronig relations, therefore we find the peak 
structure-like features as zero-crossings in the real part. 
The apparent symmetry of the self-energies is due to the 
kernel we have used for the analytical continuation pro- 
cedure, starting from the imaginary part of the Green's 
function in Matsubara frequencies. Since we have 



S'(ia;„) — duj 



Im(E'(cj)) 



lUJn 



(7) 



the imaginary part is 



which is non-vanishing only on the impurity site - in 



Im (£'(»..„)) = -^n / dJ^^^^, (8) 



which is symmetric in u). The temperature dependence 
of the self-energy documents the crossover from a single 
peak to a three peak structure in the impurity spectral 
function shown in [Fig. 7| At temperature scales above 
the Hubbard (7, correlations effects are not important 
and the self-energy essentially vanishes such that the im- 
purity spectral function reduces to the non-interacting 
one with a single central peak pinned at the Fermi en- 
ergy due to particle-hole symmetry. Lowering the tem- 
perature, we observe the formation of three zero crossings 
in Re(S). Two are roughly located at a; = ±(7/2, are 
heavily damped since Im (S) is large for those frequen- 
cies, and correspond to the upper and lower Hubbard 



6 









a U = L 
□ U = 2 






. U = i 














spectral functions of Fig. 7 For U = 2we have an inverse 
Kondo temperature of /3k ~ 2.1 which is consistent since 
somewhere in the range /3 — 2 and 13 — 10 the Kondo 
resonance starts to build-up at a; = 0. For U = 4 we get 
— 50 which is again consistent with the data for the 
spectral function. 



(a) 



(b) 



FIG. 9. (Color online) The data for estimating the Kondo 
temperature Tk for the parameters V — 0.5 and A = 0.1. 
Since Xzz is a universal function with Tk as the only param- 
eter we show in (a) a data collapse of the data for different 
values of U. The data points roughly lie on the same func- 
tion, with the U = 1 points extending to the lowest relative 
temperatures, (b) shows in a log-plot the obtained Kondo 
temperatures with its dependence on U. It shows an expo- 
nential dependence on the interaction U that is expected for 
the symmetric Anderson model. 



features. The central zero crossing corresponds to the 
Kondo resonance. It has a narrow line shape and hence 
has a small value of Im (E). That these central crossings 
have the same slope is somehow expected, since the slope 
should be proportional to the Kondo temperature of the 
systerrPS. 



A. Estimating the Kondo temperature 

To get an estimate of the involved Kondo tempera- 
tures Tk we carried out a data collapse in |Fig. 9a| of the 
susceptibilities shown in |Fig. 6h\ We find the numerical 



values for the inverse Kondo temperature /3k — T, 
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from [Fig. 9a| It is known that the susceptibility of the 
Kondo model should be a universal function 



TkXzz = F 



T 
TV 



(9) 



with Tk as the only scaling parametei^^s!. With the 
Kondo temperatures at hand we performed a cross check 
of the obtained values of Tk and assumed the validity of 
the asymptotic behaviour of the Kondo temperature for 
the symmetric Anderson model, given by 



Tk oc e 



(10) 



with the density of states po^^ ■ Note that for the Kane- 
Mele model in the considered parameter range, the Fermi 
velocity of the edge state is set by the spin orbit cou- 
pling. Hence po oc 1/A. The straight line in the log-plot 
of |Fig. 9b| indeed confirms this behaviour. Additionally, 
another cross-check is available by means of the impurity 



B. High temperature regime 

The high-temperature regime is defined by the lack 
of any visible structure in the dot's spectral function. 
/? = 0.1 in Fig. 7a is a good example, it shows just some 
lorentzian peak around w = 0. With the presence of the 
small parameter f3U it is obvious that any interaction 
induced correlation effects to the spectral functions are 
negligible. Therefore in this regime the spectral func- 
tion is that of the uncorrelated system. Since all corre- 
lation effects are thermally washed out, the notion of a 
self-energy is meaningless and hence the self-energy con- 
tribution C„(r, a;,E(a;)) in Eq. (3) vanishes. Therefore, 



the lattice spectral functions look indistinguishable to the 
non-interacting case |Fig. 5a[ Since we are considering 
the particle-hole symmetric point, the occupancy of the 
impurity site is pinned to half-filling. Hence in the ab- 
sence of interactions, the double occupancy (n'^nf) takes 
the value 0.25. As apparent in |Fig. 6a| this value is ap- 
proached as /? — 0. Finally, in this high temperature 
limit the spin susceptibility shows a l/T behaviour as 
apparent from |Fig. 6b[ 



C. Local Moment regime 

The formation of a local moment is at best character- 
ized by the quantity 



{n1}{nf} 



(11) 



Since as mentioned previously the denominator of the 
above equation is pinned to 0.25 by particle-hole symme- 
try, the formation of the local moment boils down to the 
suppression of the double occupancy. We will associate a 
characteristic energy scale of this regime by the dip in the 
double occupancy of |Fig. 6a[ The local moment regime 
corresponds to the regime where the Hubbard bands at 
develop in the dot's spectral function. In Fig. 7a 
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/3 = 2 is a good example of that. At these intermediate 
temperatures we have one occupied spin state below the 
fermi energy at — ^ with a single electron. This gives rise 
to essentially a free spin- i degree of freedom: a local mo- 
ment. The energy scale at which the double occupancy 
is enhanced before saturating marks the onset of the su- 
per exchange scale. This scale is set by /U for our 
particle-hole symmetric impurity problem. In the con- 
text of the helical liquid the spin flip scattering generated 
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FIG. 10. (Color online) The spectra for the local moment regime a,t j3 = 2,U = 2,V = 0.5 and A = 0.1. (a) shows all spectral 
functions around the impurity, whereas (b) shows a cut of the spectral functions along r = Q. (c) and (d) show the self-energy 
contribution C„(r, tj, E(tj)). C„(r = Q,uJ, is shown in (d). A neat thing is that the contribution C due to the self-energy 

seems to exactly cancel the effects of B, the hybridization of the impurity with the bath, since the contribution in e.g. |Fig. lOd] 
has opposite sign of |Fig. 5d| 



by the super-exchange scale corresponds to single parti- 
cle back scattering. These processes will hence reduce 
the conductance as noted in ref.^iSl. In the local moment 
domain, the site resolved single particle spectral function 
of the edge and bulk states in [Fig. 10a] show no sign of 
Kondo screening. In particular there is no deflection of 
the edge current perceivable. |Fig. 12b| shows that the 
change of v4„(r, a;) relative to some very distant point of 



reference has a very small amplitude. Due to Eg. (3)| this 
means an approximate cancellation of the hybridisation 



effect from |Fig. 5d with the now present self-energy ef- 



fect shown in Fig. 10d[ The edge state is in that sense re- 
stored to the true boundary of the lattice as if no impurity 
were present. This confirms the robustness of the edge 
state to a free local moment at least in view of the spec- 
tral functions. Nevertheless we expect the conductance 
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through this edge state to decrease in this regime due 
to the possible backscattering spin-flip processes which 
is now available due to the impuritj^^. Lowering the 
temperature further should break this match and deform 
the edge state again, with the emergence of the Kondo 
resonance at cj = in the impurities spectral function. 



D. Kondo regime 

Lowering the temperature further we enter the Kondo 
regime. The prominent signature of this regime is the 
emergence of the Kondo resonance in the dot's spectral 
functions as seen for example in |Fig. 7a| at /3 = 200. The 
Kondo resonance shows up as a dip in the local spectral 
function of the edge state electrons at the position of the 
impurity {r — and n = 0) as seen in Fig. lib The ori- 
gin of these spectral features lies in the formation of the 
Kondo singlet which entangles the impurity spin with the 
spins of the surrounding electrons. On an energy scale 
set by the Kondo temperature we expect the impurity 
site to act as a potential scatterer and hence lead to the 
deflection of the edge current into the bulk. As apparent, 
the dip in the spectral weight at rt = in Fig. lld| [b) is 
accompanied by the emergence of a peak at n = 1 and 
n = 2. At this very low temperature scale, T/Tk ~ 0.01, 
it is this deflection of the edge state discussed by Ma- 
ciejko et al.^'^ which restores the conductance to unitar- 
ity in the limit T ^ 0. Although we have decreased the 
temperature by two magnitudes the form of the change 
in |Fig. lid] still looks the same, only its magnitude has 
decreased. But of course now the self-energy term C and 
the hybridization term B do not match anymore as nicely 
as for /3 = 2. Hence there is not a full cancellation of the 
hybridization effect. A further decrease of the tempera- 
ture should bring us deeper into the Kondo regime with 
a deflection of the edge state around the impurity as pre- 
dicted by Maciejko et al.l^. This confirms that, to the 
outside world, the singlet state of a magnetic impurity 
and bath electrons has the same low-energy features as 
plain potential scattering. We only expect this circum- 
vention of the edge state on an energy scale set by the 
Kondo temperature Tk, since beyond this energy scale 
the equivalence to a potential scatterer is not tenabl^^H, 
Note that since the width of the Kondo resonance is of 
the order of TK^ the dip in the bulk spectral functions is 
of the same size. On the other hand, the width of the 
edge state is given by the spin orbit coupling A > Tk- 
|Fig. 12| shows the change attributable to the impurity 
measured against a very distant point in the same orbital 
n: A„(r, w) = ^„(r, w) — A^{ijj). We see the strong local 
effect of a potential scatterer in [Fig. 12a| In [Fig. l"2b| we 
see that a correlated impurity aX P — 2 has a negligible 
effect on the bath. This corresponds to the local moment 
regime. As the temperature drops below the Kondo scale 
(s ee the dat a sets at /3 = 7 in Fig. 12c and &% j3 — 200 
in 



VI. SPATIALLY RESOLVED DOT BULK SPIN 
SPIN CORRELATION FUNCTIONS 

To provide a different point of view on our study of 
the Kondo cloud in the bath system that does not rely 
on an analytic continuation procedure we now turn our 
attention towards the site resolved spin spin correlation 
functions 



between the impurity spin — ^ 
located at a particular site S'^ (r, n 



(12) 



and a spin 



This enables us to define the Kondo cloud as the region 
of substantial entanglement of the impurity spin with 
a particular bath site. This is complementary to our 
previous results where we have defined the Kondo cloud 
as the region where the edge state is suppressed. 



A. A 2D overview 

First, we consider the non- interacting case and find 
non-negligible correlations that are confined to the edge 
of the system. Along the edge the spin spin correlations 
decay as as already pointed out in Ref.^^ for a ID 
system of electrons. Clearly, this power law holds only 
in the zero temperature limit and at finite temperature 
an exponential decay sets in beyond the thermal length 
scale c>c This similarity to the one- dimensional 

case provides yet another confirmation of the ID nature 
of the edge state. 

The local moment regime just shows a small amount of 
correlation in |Fig. 15| since due to the high-temperature 
all long-range effects are destroyed. This is in contrast 
to the behaviour of the spectral functions in [Fig. 10a| 
where we see no signature of the impurity in the bath 
system. If we now lower the temperature into the Kondo 
regime in 



Fig. 13 we see that the effect of the impurity 



Fig. 12d) the same defiection of the edge current as 



observed for the potential scatterer emerges. 



mostly extends into the helical liquid in the lower edge 
and develops some spatial structure. 



B. Correlation functions along the edge 

We now focus our attention on the correlations along 
the edge as a function of temperature, Hubbard inter- 
action and spin-orbit coupling A. Borda et al.'^ have 
studied the spatial behaviour of the spin spin correlation 
functions of an Anderson impurity embedded in a one- 
dimensional wire as bath system. They observed at a dis- 
tance of « vfPk {f^F is the Fermi velocity) a crossover 
from an r^^ behaviour to an behaviour. For finite 
temperatures they predict the onset of an exponential 
decay at ~ vfP- Their study shows that the spa- 
tial decay is oscillating with a wave vector k (x kp — n. 
It can already be guessed from the 2D overviews that 
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FIG. 11. (Color online) For the parameters /3 — 200, U — 2,V = 0.5 and A — 0.1, (a) shows a frontal view on the full spectral 
functions. The dip at the impurity is visible as well as the progression to the full edge state at the edge, (b) shows a cut of the 
same spectral functions now for increasing n. The displacement of the edge into the bulk is visible, (c) and (d) show the effect 
due to the self-energy Cn(r, w, S(aj)). Especially figure (d) shows that in comparison with Fig. lOd the functional form of C 
seems to be the same although we have increased /3 by two magnitudes, but the amplitude is reduced from about 0.6 to 0.4. 



our system does not show oscillations, which is consis- 
tent since cos{2kp) = 1. Fig. 14 shows that the general 



trend of these predictions made for a ID chain of elec- 
trons also holds for the ID helical liquid of the edge state 
of a topological insulator. In Fig. 14 (a) we find the de- 



pendence on the Fermi velocity, which is expected to be 
proportional to the spin-orbit coupling A. In the plot for 
A = 0.2 the expontial falloff is shifted outside of the visi- 



ble part of the edge channel, although the plot for A = 0.1 
and A = 0.2 are at the same temperature. In this plot we 
already have introduced a couple of lines that are meant 
as guides to the eye. The yellow dashed lines denote the 
exponential decay that sets in beyond the thermal length 
scale £,T oc VfP. The straight dashed magenta lines de- 
note the power law decay that is present at distances 
£,K > r > S^T- Finally, the straight dashed orange lines 
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FIG. 12. Here we show the changes to the spectral functions 
in the vicinity of the impurity relative to some very distant 
reference point that feels no effect of the impurity. This cor- 
responds to the quantity A„(r, a;) = A„{r,u)) — A°(tj). (a) 
shows the changes due to an uncorrelated impurity, (b) is 
the change in the local moment regime at /3 = 2. Note the 
small amplitude of A, which means that in this regime the 
bath feels only a negligible effect of the impurity, (c) is in 
the Kondo regime for 13 — 7. (d) is in the Kondo regime for 
/3 = 200. Note that only (b) - (d) share the same scale on the 
A„(r, uu) axis. 



denote the decay that is present for r < S^^- We see 
that for A = 0.2 the crossover from a r^^ behaviour to a 
decay is approximately shifted from r w 7 to r w 16 
and the thermally induced exponential suppression of the 
spin spin correlation happens much later. We can esti- 
mate the thermal cutoff scale by fitting exponentials (the 
yellow dashed lines), e~'^ , to the tails of the plots for 
/3 — 100 and we find a consistent value of ~ 8.7, for all 
values of U. A more detailed analysis of the temperature 
dependence for the point [/ = 2 is found in [Fig. 15[ In 
Fig. 14| (b) and (c) we can compare the temperature ef- 
fects for [/ = 3 and we see that at twice the temperature 
the exponential decay is not visible anymore. Compar- 
ing the plots from (a) to (d) we can trace the shift of the 
cross-over from an behaviour to an r^^ decay with 
increasing U and therefore with the Kondo temperature. 
Fig. 14 (c) also shows the independence of the cross-over 
point of the algebraic decay due to the Kondo effect with 
respect to the external temperature. 
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FIG. 13. (Color online) A 2D overview of the color-coded 
correlation function {S^S^ {r,n)) with the full spatial depen- 
dence at /3 = 100, A = 0.1 and U = 2. The impurity is the 
green spot at (r, n) — (0, 0) in the middle of the lower edge. 
Positive values are shades of green, negative values are shades 
of red and neutral values are blue. The color coding is linear. 
Note that this does not mean a linear perception of the color 
values. The comparatively long chain of sites having various 
shades of red at the bottom of the diagram is the extent that 
the correlation reaches into the edge state. The correlation 
extends comparatively far into the edge state but is almost 
immediately suppressed away from the edge state. 



VII. SUMMARY 

We have studied a magnetic impurity coupled to a he- 
lical edge state as modeled by a Kane-Mele Hamiltonian 
on a slab geometry. Due to time reversal symmetry the 

has 



effective action of the impurity orbital (see Eq. (5) 



precisely the same form as the generic SIAM such that 
the local physics is identical. In particular the Hubbard 
scale marks the appearance of a local moment which cou- 
ples magnetically via the superexchange scale J oc V^/U 
to the conduction electrons. Below the Kondo temper- 
ature, Tk, the magnetic moment is screened due to the 
formation of an entangled singlet state of the magnetic 
impurity and conduction electrons. We have shown these 
commonalities numerically with the double occupancy, 
the spin susceptibility and determined the Kondo tem- 
perature with a data collapse. 

The differences to generic Kondo physics are non-local. 
As shown in Ref.^^, if the interactions along the helical 
edge are not too strong, spin-flip single-particle backward 
scattering processes are expected to be irrelevant, such 
that in the low temperature limit the conductance should 
reach the unitarity limit. The mechanism which allows 
this to occur is the deflection of the edge state into the 
bulk, thus avoiding the Kondo singlet. By computing 
the temperature dependence of the site resolved density 
of states by an analytic continuation of the impurity self- 
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FIG. 14. (Color online) The spatial dependence of the corre- 
lation between the impurity's spin and the spin of a site of the 
edge \{S^S''{r, n))\ with a hybridization ofV = 0.5. Diagram 
(a) shows the dependence on the Fermi velocity at (7 = 1 of 
the cross-over point from an to an behaviour, ^x, as 
well as of the point where the decay crosses over into an 
exponential law. Note that for = vfPk both the Fermi 
velocity and Kondo temperature depend on A: vf oi \ and 
pK oc exp ( j) . In the graph for A = 0.2 the exponential falloff 
is shifted outside of the visible lattice although the temper- 
ature is kept constant. The dashed lines are guides to the 
eye and are explained in the main text. Diagrams (b) - (d) 
are at A = 0.1. In diagram (b) we show data for higher Tk 
aX U = 2. We see that for the higher temperature P = 200 
the thermal decay at the end of the plot is shifted outside of 
the visible part of the lattice and we can clearly identify the 
regimes with power law like behaviour. In the diagram for 
P — 100 the decay is already dominated by the thermal 
decay. Figure (c) and (d) show that further increments of U , 
which gives higher values of /3x, shifts the cross - over point 
to larger distances. 



(S3S'(r,„)) 

(a) 




(b) 



FIG. 15. (Color online) (a) shows the spin spin correlation 
functions in the local moment regime. Due to the rather high 
temperature we see that the correlation effects are quickly 
suppressed. Additionally, the absolute value of the correla- 
tions is an order of magnitude lower. The right figure pro- 
vides a more detailed view of |{5'^S^(r, 0))| aX U = 2, but at 
different values of /3 in a logarithmic plot. Here we can notice 
the mirror symmetry around the impurity. We see that for 
P = 10, which is in the l ocal moment regime (compare the 
spectral functions Fig. 7a I, the decay of the correlation func- 



tion is immediately exponential. This is consistent with the 
quick suppression in the left figure. Decreasing the tempera- 
ture to P = 100 we see that the exponential decay sets in at 
around ~ 8 whereas for P — 200 at ~ 16. 



spatial dependence of the spin-spin correlation functions 
between the local moment and conduction electrons. At 
low temperatures the spin-spin correlations do not extend 
significantly into the bulk, and exhibit a power law decay 
along the edge. In particular, as a function of tempera- 
ture, we can observe the thermal cutoff scale oc vjP 
beyond which exponential decay sets in, as well as the 
characteristic cross-over scale around oc Vf/Tx from 
an to an law. This cross-over scale provides a 
measure of the Kondo screening cloud. Of significant 
interest is an explicit calculation of the temperature de- 
pendence of the conductance along the edge. In partic- 
ular, in the local moment regime, where we observe no 
deflection of the edge state, one expects a decrease of the 
conductance due to back-scattering spin flip processes off 
the impurity spin. Below the Kondo scale, the deflection 
of the edge state along the new boundary of the system 
- as defined by the topology of the Kondo cloud - should 
restore the conductance to it's unitarity limit. 
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Appendix A: Analytical continuation 

The analytic continuation of the Matsubara Green's 
function G(n, fc,iu;„) to the real frequency axis is a no- 
toriously hard problem and requires potentially large 
amounts of computer time. For Monte Carlo data, 
experience shows that the most reliable spectra can 
be obtained using the stochastic maximum entropy 
method^i'^ for analytic continuation. This method, 
however, makes use of a Monte Carlo simulation to find 
the best spectrum and has to be performed for every 
pair of indices {n, k) individually. Clearly, this increases 
the computational effort substantially. In order to cir- 
cumvent this problem, we propose to adopt an idea that 
has been sucessfuUy used in the context of the dynamical 
cluster approximation (DCAJP^l Instead of the Green's 
function, we calculate the analytic continuation of the 
self-energy. The reason why this is beneficial is that the 
self-energy matrix has only one non-vanishing and diag- 
onal 2x2 spin block S(z) corresponding to the impurity 
site and leads to the special form of Dyson's equation: 



^G%^,{z) G'ljz) 
yGlA^) GUz 



G^z) GUz) 




Gr,r{z) Gp,d{z) 

GdA^) Gd^iz) 



(Al) 



Obviously, the knowledge of S](z) on the real axis is there- 
fore sufficient to obtain the full interacting Green's func- 
tion matrix G{ijJ + irj) for the whole lattice as the nonin- 
teracting Green's function can be calculated exactly 
with little difficulty by virtue of the resolvent formalism. 
Solving equation | ( A 1 ) | for Yi{z) yields 



As Gd,d{}'-^n) can be calculated in CT-INT, we can there- 
fore obtain E(ia;„). In order to analytically continue 
to the real axis, we have to study its asymptotic be- 
haviour for large frequencies. Starting from the asymp- 
totic series for G^''^(icj„) 



oo (0) 

— a). 



we obtain through inversion: 



(A3) 



E(ia;„) = (02 - 0°) + 



1 CL\ — {a^) -1-0103 — 0103 
On a I 



O 



(A4) 



This result can be obtained by truncating equation |(A3)| 
at different orders and one indeed finds out that higher 



terms of the Green's function's asymptotic serie^^do not 
contribute to the first two terms of the self-energy. 

In order to employ the stochastic maximum entropy 
method for 'S{z) directly, we introduce a slightly different 
quantity as already shown in reference 1331 



E'(z) 



[l](z) - (02 - O^)] Oi 



(A5) 



This quantity has exactly the same properties as the 
Green's function itself, namely that its asymptotic se- 
ries starts with the corresponding spectral function 

has a sum rule J da;v4x;(a;) — tt and that it does not have 
a constant term. 

In principle, these properties could be corrected for in 
the maximum entropy procedure but the quantities 02 
and 03 can only be obtained up to a statistical error- 
bar and therefore the correct inclusion of these errors is 
very cumbersome. Performing the transformation [(A5) | is 
therefore a very straightforward procedure, as the thor- 
oughly bootstrapped covariance matrix of Yl will contain 
all uncertainties stemming from the CT-INT calculation. 

The calculation of the constants in the asymptotic se- 
ries of the self-energy is a straightforward calculation of 
moments^^of the spectral function A(u}) corresponding to 
the Green's function Gd,d and yields for Ag, (w) 



U 
2^' 



02-0°=C/(d^_,d_,) 

(0^)^-0^-1-01(03-0!]) _ 

Ol 

At half filling, (dj^do-) = |, so the constant term 02 



(A6) 



(A7) 



(A2) of the self-energy vanishes. 



Appendix B: Diagonal Gd,d due to TRI 

Since [iJ, T] = 0, the operators H and T have a com- 
mon basis in which they are diagonal. In the case that 
T'^ — —1 holds (which is the case if the total spin is of 
the half-integer type) we have for every state el = | o) an 
orthogonal state 62 = ^jo). In these basis states T has 
the matrix representation 




(Bl) 



where 1 denotes the identity matrix acting in the respec- 
tive sub-sector of the state. By that notation we have 
essentially just renamed the states of the Hilbert space. 
T is diagonalized by 



U = 



1 



1 il 



V2 \il 1 



(B2) 
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with two eigenvectors named | + ) and | — ) . U also block- 
diagonalizes H and hence the Green's function matrix is 
block-diagonal in the |+) and | — ) basis since propagation 
with H will not generate matrix elements between the 
othogonal sub-spaces. 



G 



G+^ 







G- 



(B3) 



Transforming G back to the original basis we have 



G = UGU"^ 




iG++ - iG- 
G++ + G- 



(B4) 



To show that the ofF-diagonals of this matrix vanish in 
the sub-space of this matrix where the impurity lives, we 
consider the impurity Green's function Gd,d in the | -I- ) 



and I — ) basis. Then, using time-reversal symmetry. 



THT-^ = H 

we have for t > 0, 

GrdM = Tr(e-^^4(r)4) 

= Tr{e-^^"^-'dl{r)d, 
= Tr (Te-^^T-ie-^^4e-^^4) 
= Tr {e''^"e-^"T-^dlTe'-"T'^d,T) 



(B5) 



(B6) 



using T ^dgT = sd-s (the impurity lacks a momentum 
quantum number) we have 



G 



d,d 



s^Tr 



{e-^"dUr)d-. 



(B7) 



Therefore we have G^'J^ = G^ thereby ensuring that 



the off-diagonals in Eq. (B4) vanish. This in turn gives 



the diagonality of an impurity Green's function Gd.d 
just due to time reversal symmetry. 
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